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Abstract. European VLBI Network spectral imaging of the "compact double" radio source 2050+364 in the UHF band at 
1049 MHz has resolved the Hi absorbing region, and has shown a faint continuum component to the North (N), in addition to 
the well-known East- West double (E, W). 

Re-examination of VLBI continuum images at multiple frequencies suggests that 2050+364 may well be a one-sided core-jet 
source, which appears as a double over a limited frequency range. One of the dominant features, W, would then be the innermost 
visible portion of the jet, and could be at or adjacent to the canonical radio core. The other, E, is probably related to shocks at a 
sudden bend of the jet, towards extended steep-spectrum region N. 

A remarkably deep and narrow H i absorption line component extends over the entire projected extent of 2050+364. It coincides 
in velocity with the [O m] optical doublet lines to within 10 km s . This H i absorption could arise in the atomic cores of NLR 
clouds, and the motion in the NLR is then remarkably coherent both along the line-of-sight and across a projected distance of 
> 300 pc on the plane of the sky. 

Broader, shallower H i absorption at lower velocities covers only the plausible core area W. This absorption could be due to gas 
which is either being entrained by the inner jet or is flowing out from the accretion region; it could be related to the BLR. 
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1. Introduction 

2050+364 (J2052+3635, DA 529; J2000 coordinates 
20 h 52 m 52.0549 s +36°35'35.300", Beasley et al. 2002), is 
one the original members of the group of radio sources 
described as "Compact Doubles" (CDs) by Phillips & Mutel 
(1981, 1982), and then studied in more detail by Mutel, 
Hodges, & Phillips (1985). Phillips & Mutel (1981) suggested 
that, in contrast to most of the radio-loud AGN imaged 
with VLBI at GHz frequencies, CDs are "mini-lobes", on 
(sub)galactic scales, at the ends of bi-directional outflows from 
an unseen central ("core") component. The spectrum of the 
two components in compact doubles is typically not flat but 
peaked at a frequency around a GHz; such sources are termed 
Gigahertz Peaked Spectrum (GPS) sources (see the review by 
O'Dea 1998). The spectral shape can often be reasonably well 
understood from synchrotron self-absorption in ~ 1 GHz radio 
emitting regions of ~ 100 mas cross-section, but there are 
indications that in some sources, the low-frequency spectral 
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shape is also determined (in part ?) by free-free absorption 
(e.g., Bicknell, Dopita, & O'Dea 1997; Kameno et al. 2003, 
hereafter K03; Risaliti, Woltjer, & Salvati 2003). 

The CD and GPS source classes are related to the group 
of Compact Symmetric Objects (CSOs), first so termed by 
Conway et al. (1994) and Wilkinson et al. (1994). In CSOs, 
a central compact core component is visible in between two 
lobe-like components, albeit usually only faintly and often 
only at high observing frequencies. Thus, in CSOs the two- 
sidedness of the radio emission is proven beyond reasonable 
doubt. Models have been considered in which the lobes are 
kept to subgalactic dimensions by a dense confining medium 
(e.g., De Young 1993; Carvalho 1994, 1998). However, in a 
number of CSOs (e.g. Owsianik & Conway 1998; Polatidis & 
Conway 2003; Gugliucci et al. 2005) lobe advance velocities 
of a few tenths c have now been observed, and this provides 
powerful evidence that CSOs structures signify young radio 
sources, some with kinematic ages as low as a few hundred 
years. 

GPS radio sources at moderate redshifts (z < 1) are typi- 
cally identified with passively evolving elliptical galaxies (e.g., 
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Snellen et al. 1998). The lifecycle of these sources is a subject 
of current interest, and is not well understood: at least some 
of them probably evolve to become fully fledged double-lobed 
FR I or perhaps even FR II radio sources, which typically have 
linear sizes of 100 kpc or even up to 1 Mpc, although, based 
on number counts, they should decrease in luminosity as they 
grow (e.g., Fanti et al. 1995; Readhead et al. 1996; O'Dea 
& Baum 1997; Alexander 2000; Snellen et al. 2000). On the 
other hand, perhaps some compact sources may have multiple 
very short-lived active episodes (e.g., Gugliucci et al. 2005). 
Compact double radio sources are also particularly promising 
objects to study in the context of AGN fuelling and radio source 
evolution models, and the interaction between the interstellar 
medium and jets, because they lie entirely within the inner 
galaxy, on kpc or even sub-kpc scales. 

One attractive possibility is to observe the 21cm hyper- 
fine line of atomic hydrogen in absorption against bright radio 
structure, to study the kinematics and density distribution of 
the atomic gas. VLBI offers spatial resolution on the scales of 
the NLR and sometimes even the BLR, not accessible by any 
other means in galaxies at any appreciable redshift. However, 
for most CPS/CSO/CSS radio sources the frequency of this H i 
line is redshifted outside the traditional observing band. But 
profiting from the advent of UHF receivers on interferometers, 
first at the Westerbork Synthesis Radio Telescope (WSRT), 
and then at many telescopes in the European VLBI Network 
(EVN), it has become possible to study the presence and dis- 
tribution of associated H i in absorption against many compact 
radio sources. 

A WSRT survey of a sample of compact radio sources 
was published by Vermeulen et al. (2003, hereafter V03), and 
further analysed by Pihlstrom, Conway, & Vermeulen (2003). 
About one third of the objects surveyed had detectable H i ab- 
sorption, and the peak depth was found to be anti-correlated 
with linear size. A 16 % peak depth HI absorption line at 
z = 0.3547 1 was discovered towards 2050+364. This is uncom- 
monly deep; for sources of similar overall linear extent, absorp- 
tion depths of 0. 1-1 % (or upper limits at that level) were found 
to be rather more typical by V03. 2050+364 was therefore in- 
cluded in one of the first sessions with the new UHF receivers 
on the EVN, in order to resolve the HI absorption region spa- 
tially. We describe the data analysis procedure in Sec. 2, and 
present the results in Sec. 3. 

2050+364 has received comparatively little attention in its 
own right during the last two decades. Perhaps this was because 
of its low Galactic latitude, b - -5°, in the Cygnus superbubble 
region (/ = 79°), which is a handicap for obtaining a secure op- 
tical identification because of crowding, and hampers low fre- 
quency radio studies of the object because its image is broad- 
ened by interstellar scattering (Dennison et al. 1984, Mutel & 
Hodges 1986). Biretta, Schneider, & Gunn (1985) and O'Dea, 
Baum, & Morris (1990) have both found a good match between 
the radio position and an m r = 21.1, (r-i) = 0.1 galaxy. A spec- 
trum was published by de Vries et al. (2000) showing promi- 

1 H = 70 kms" 1 Mpc" 1 , C2 m = 0.27, and Q A = 0.73 are adopted 
throughout this paper, so that an angular size of 1 mas corresponds to 
5.0 pc at z = 0.3547. 



nent H/3 and [Oin] emission lines. The latter have a redshift 
z = 0.3547; we will critically review the redshift of 2050+364 
in Sec. 5.2 and the optical identification in Sec. 5.3. Meanwhile, 
the low galactic latitude of 2050+364 has been turned into an 
advantage by using it as a background probe for interstellar 
scattering (e.g., Fey & Mutel 1993). This has culminated in a 
multi-frequency VLBI imaging study by Lazio & Fey (2001, 
hereafter LF). Another multi-frequency VLBI dataset is avail- 
able from the investigation into possible free-free absorption in 
a sample of GPS sources by K03. Their results are described in 
Sec. 4, and used in our discussion on the nature of 2050+364 
in Sec. 5. 

2. Observations and data processing 

2.1. Observations 

On 1999 September 08-09, in one of the earliest sessions at 
UHF frequecies with the European VLBI Network (EVN), 
14 hours were spent observing 2050+364, and the cali- 
brators DA 406, 3C 454.3, and 3C84. Four telescopes pro- 
duced usable data: dual circular polarisations were available 
at Effelsberg and the WSRT, while Jodrell Bank had only 
LCP, and Onsala recorded dual linear polarisations. The (u,v)- 
coverage on 2050+364 is shown in Fig. 1. The observing band, 
1046.6 MHz to 1050.6 MHz, was centred on the absorption 
line discovered with the WSRT by V03. A more extensive de- 
scription of early EVN UHF observing characteristics and data 
analysis procedures will be given in Vermeulen et al. (in prepa- 
ration). 




-2 
U (10 6 A) 



Fig. 1. The (u,v)-coverage obtained on 2050+364 with the EVN tele- 
scopes Effelsberg, Jodrell Bank, Onsala, and WSRT. Crossing points 
provide important self-calibration constraints. 
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2.2. Initial calibration and data averaging 

The initial data processing (fringe fitting, spectral passband cal- 
ibration, and a priori complex gain calibration) took place in 
the NRAO Astronomical Image Processing Software (AIPS), 
using the calibrator sources. After time-interpolated transfer of 
the solutions, and residual fringe-fitting on 2050+364, the data 
were averaged into 256 independent spectral channels (each 
4.5 kms~' wide), and 60 sec time samples. This increased the 
sensitivity per visibility. The target source, 2050+364, was de- 
tected on all baselines. 

All available polarisation products were also averaged to- 
gether, including the cross-correlations between the linear po- 
larisations from Onsala and the circular polarisations from the 
other telescopes. We believe the sensitivity gained to total in- 
tensity is more important than any possible resultant limitations 
on absolute flux calibration accuracy or on (image or spectral) 
dynamic range; we think these are more affected by the sparse- 
ness of the array and the lack of complete system temperature 
data and gain curves. 




Right Ascension (mas) 



Fig. 2. The dirty beam obtained with the (u,v)- 
coverage shown in Fig. 1; contour levels are drawn at 
-0.2, -0.1, -0.05, 0.05, 0.1, 0.2, 0.4, and 0.8. This beam is 
benign for deconvolution of the continuum structures in Fig. 3 

Much care was taken to obtain consistent visibility calibra- 
tion, for the rather sparse array, and in the occasional presence 
of external radio interference (although its impact is mitigated 
since it is typically not coherent over VLBI baselines). Relative 
antenna gains were derived from requiring consistent visibil- 
ity amplitudes on the calibrator sources, while the overall flux 
density scale was set such that in the final self-calibrated im- 
age (see below) the total flux density is equal to 3.34 Jy at 
1049 MHz, a value derived by interpolation from total flux den- 
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sities at a number of other frequencies (Salgado et al. 1999; 
White & Becker 1992). We estimate the overall uncertainty in 
the flux scale of this UHF VLBI dataset to be as much as 20 %, 
but the main astrophysical results are unaffected, since they de- 
pend on opacities rather than on absolute flux densities. 
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Fig. 3. Continuum VLBI image of 2050+364 at 1049 MHz. The 
restoring beam is 50 x 38 mas in position angle -27°. The peak flux 
density is 2.42 ly/beam; contours start at 0.015 Jy/beam (3 times the 
r.m.s. fluctuation level away from the continuum structure), and in- 
crease in factors of 2. The absolute flux density scale is uncertain by 
20 %, but this does not affect relative flux densities; see Sec. 2.2. The 
three circular Gaussian components used for modeling (Table 1) are 
overplotted. 



2.3. Self -calibration; continuum imaging and modeling 

The Caltech DIFMAP software package (Shepherd, Pearson & 
Taylor 1977) was used for all further calibration and analy- 
sis. A continuum dataset was formed by averaging all of the 
line-free channels on both sides of the H i absorption line, of 
which the location and width were taken from the WSRT sur- 
vey (V03). The available VLBI array was sparse, but contained 
some crossing points (see Fig. 1), and careful self-calibration 
was feasible. A double source structure was already evident in 
the first image made. After initial phase self-calibration, some 
additional flux density to the north-east quickly became appar- 
ent. At that time, we were unaware of the images published 
by LF, which at the lower frequencies also show emission in 
that area. Extensive tests of the reality of this third component 
were convincing: our data demand its presence. There are no 
major sidelobes of the dirty beam (shown in Fig. 2) to affect its 
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Table 1. The three circular Gaussian components used for modeling. 
Their continuum flux densities are listed. The absolute flux density 
scale is uncertain by 20 %, but this does not affect relative flux densi- 
ties; see Sec. 2.2. These components are overplotted on the image in 
Fig. 3. 



Component 


Flux 


Radial 


P.A. 


Diam 




(Jy) 


dist (mas) 


(deg) 


(mas) 


W 


0.79 








13.1 


E 


2.44 


57 


99 


5.2 


N 


0.11 


74 


66 


7.7 



deconvolution with respect to the brighter parts of the source. 
We have performed careful iterative cycles of self-calibration 
with imaging/cleaning/modeling, both in trial runs where we 
kept excluding the third component from the models, and in 
runs where we admitted the presence of this additional contin- 
uum feature. The latter both gave cleaner-looking images and, 
more significantly, consistently required less extreme and more 
time-stable self-calibration coefficients. 

In order to restrict the number of free parameters, the 
sky model fitted to the visibility data during iterative self- 
calibration cycles consisted of three circular Gaussian compo- 
nents: W(est), E(ast), and N(orth); their parameters, fitted to 
the visibility data, are given in Table 1 . The final self-calibrated 
continuum visibility data were used to produce the cleaned and 
restored image displayed in Fig. 3; symbols showing the three 
model components are overplotted. 

2.4. Spatially resolved spectroscopy 

The sensitivity of the data and the sparseness of the array did 
not allow generation of a full spectral image cube of accept- 
able signal-to-noise ratio. Instead, spectra of the W, E, and N 
areas of the source were determined by re-fitting the flux den- 
sity of each of the three Gaussian model components derived 
for the line-free continuum, but then separately for each of the 
256 spectral channels. The positions and diameters of the three 
components were not varied; they were fixed at the values de- 
rived from the continuum. 

The resultant spectra are shown in Figures 4 and 5, as opac- 
ity with respect to the continuum strength of the appropriate 
components; this avoids the uncertainties in the absolute flux 
density scale, discussed in Sec. 2.2. Zero velocity for Figs. 4 
and 5 corresponds to the redshift of [Oiii], z = 0.3547, pub- 
lished by de Vries et al. (2000); this choice will be further dis- 
cussed in Sec. 5.2. Note that a less precise value of the optical 
redshift was used by V03. 

Gaussian line profiles were fitted to the most obvious fea- 
tures of the Hi spectra at the three locations, using Origin 6.1. 
The fitted profiles are overplotted on the data in Fig. 6, and their 
parameters are listed in Table 2. These serve only to obtain a 
rough quantification of the absorber properties; in particular, no 
attempt was made to find an optimally fitting set of Gaussians 
to cover the complex observed line spectrum at W. 
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Fig. 4. The H i absorption spectra, expressed in opacity, for the W, E, 
and N radio components of 2050+364. Zero velocity is chosen to be 
the [Om] redshift, z = 0.3547 (see Sec. 5.2). 
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Fig. 5. Enlargement of part of Fig. 4, to allow detailed comparison of 
the H i opacity profiles towards W and E. 

3. Results 

3.1. VLBI continuum structure of 2050+364 at 1 GHz 

The 1049 MHz continuum VLBI image shows the W and E 
components, which constitute the well-known "classical dou- 
ble". At this frequency the W:E flux density ratio is close to 1:3; 
this is constrained better than the flux densities themselves. The 
two components are separated by about 57 mas (285 pc) along 
a PA of ~100°, which matches the relative locations found at 
other frequencies by LF. The component sizes we find at 1 GHz 
are only slightly larger than those found at 1 .67 GHz by LF; we 
believe that the accuracy of our fitted component sizes is lim- 
ited by the relatively poor resolution of our array (synthesised 
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Fig. 6. Fitted Gaussian line profiles overplotted onto the H i absorption 
spectra of the three regions. The parameters are listed in Table 2. 



beam 50 x 38 mas), although it is interesting that the angular 
diameter we find for E, 5.2 mas, is comparable to the scattering 
disk diameter determined by LF, 6-8 mas (as discussed below, 
this could be somewhat of an over-estimate since LF did not 



Table 2. The parameters of the Gaussian line profiles fitted to the ob- 
served H i spectra of the three regions, as shown in Fig. 6. The derived 
column depths assume uniform coverage and T s „ = 100 K. 



Region 


Centre 


FWHM 


Peak opacity 


L * ri I 




(kms" 1 ) 


(kms" 1 ) 




10 20 cm- 2 


w 


-11.7 ±0.1 


10.3 ± 0.3 


0.212 ±0.010 


4.2 ± 0.2 




-35 ± 0.7 


5.8 ±2.1 


0.037 ± 0.022 


0.4 ± 0.3 




-71.0 ±2.3 


55.0 ± 6.9 


0.036 ± 0.006 


3.9 ± 0.8 




-128.8 ±3.6 


20.7 ± 9.2 


0.013 ± 0.008 


0.5 ± 0.4 




-166.1 ± 1.9 


24.1 ±4.8 


0.026 ± 0.007 


1.2 ± 0.4 


E 


-5.1 ±0.1 


13.4 ± 0.2 


0.263 ± 0.004 


6.8 ±0.1 




-42 .4 ± 0.9 


41.1 ± 2.4 


0.034 ± 0.003 


2.7 ± 0.3 


N 


-14.2 ± 5.1 


51.5 ±12.1 


0.091 ± 0.028 


9.1 ±3.6 



recognise sub-components which exist at the lowest frequen- 
cies. 

We have also discovered with our data (before finding it 
in the low-frequency images published by LF) an additional 
emission region to the north/northeast, which we call N. While 
we are confident about the existence and approximate location 
of N, about 74 mas (370 pc) from W in PA ~ 66°, its flux 
density (~0.1 Jy) is subject to considerable uncertainty, both in 
view of the overall flux scale uncertainty and because, with the 
sparse imaging array, some extended emission may have been 
missed. The relatively small array of EVN antennas used did 
not allow the presence of substructure within any of the three 
components to be investigated. 

3.2. H i absorption distribution 

The VLBI data reveal that the two most prominent H i absorp- 
tion features (16 % and 4 % integrated opacity) found towards 
2050+364 in the WSRT survey (V03) do not cover the sub-kpc 
scale radio source uniformly. The integrated (WSRT) profile is 
dominated by the absorption towards E, because that compo- 
nent is much brighter than the other radio source components. 

The distinct, deep absorption feature at the high velocity 
end of the absorption spectrum shows a ~1 kms -1 centroid 
velocity offset between W and E; see also Fig. 5. This ob- 
served offset is highly significant, given the velocity resolution 
of 4.5 km s _1 , and the high signal-to-noise ratio of these absorp- 
tion features. The line is also narrower at W than at E. At N, 
a corresponding absorption feature is also visible in the VLBI 
data (see Figs. 5 and 6), but, due to the low background con- 
tinuum flux density, the centroid velocity and the FWHM are 
too uncertain for a useful comparison. The peak opacity at N 
appears to be lower than at either W or E. 

The VLBI data show that the second absorption feature tab- 
ulated by V03 for the integrated WSRT spectrum (where it cov- 
ers the observed velocity range -30 km s _I to -60 km s _I with 
respect to the [Oin] redshift) has a peak opacity of t = 0.035 
at both W and E. But the detailed profile differs substantially 
between W and E; it is particularly irregular at W. The inte- 
grated profile shown in V03 is dominated by E, because it is the 
brightest feature. While the absorption at N was modeled with 
just a single Gaussian line, it is wider than the main component 
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at E and W (see Table 2, and may therefore also encompass 
some secondary absorption features. 

The VLBI data also show that there are further absorp- 
tion features, extending to an observed velocity of nearly -200 
kms -1 , at comparatively low opacities (a few percent, see 
Fig. 6). These are all detected towards W only, and not towards 
E; towards N the signal-to-noise ratio is too poor to establish 
or delimit the presence of similarly shallow features. Since E 
dominates the total flux density, these features are not readily 
apparent in the integrated spectrum of V03. 

4. Multi-frequency Radio Continuum Structure 

The compact structure of 2050+364 was imaged with VLBI at 
multiple frequencies by both LF and K03. Their projects were 
aimed at studying interstellar scattering, or free-free absorption 
statistics, respectively, and the structure of 2050+364 itself was 
not analysed in detail. However, we give our own description 
of the images here, because we will then show in Sec. 5 that 
together with our own results, a surprising new interpretation 
emerges: 2050+364 could well be a one-sided core-jet source ! 

To support the discussion, we show in Fig. 7 the flux densi- 
ties at multiple frequencies as tabulated by LF and K03, and the 
flux densities at 1 GHz from our own data (Table 1). We have 
summed over all sub-components which were given the same 
label (W or E, respectively) by LF, except, as discussed below, 
at 0.61 GHz, where we recognise one of the components as N, 
and at 0.33 GHz, where we think the single flux density listed 
by LF should in fact be ascribed to a combination of all compo- 
nents. LF and K03 both have data at 2.3 GHz and 8.4 GHz. The 
morphologies in their images appear to match well, and while 
their tabulated flux densities differ at the 20 % level, the uncer- 
tainties do not obscure the trends on which we wish to focus; 
we simply use the average of the two available measurements 
in Fig. 7, with a bar to show the range between them. 

The northeastern region, which we have called N, is clearly 
distinguishable in the 1 .67 GHz image of LF, where it consists 
of rather extended emission at a distance of roughly 30-55 mas 
from region E, to the north and northeast, in position angles 
roughly between 0° and 30° as seen from E; the position angle 
of region N as seen from W is about 65°. Note that this cor- 
responds well with the position of N in our data (see Table 1 ; 
with respect to E, the distance of N at 1049 MHz is 41 mas, in 
a position angle of 16°). There are, understandably, no model 
components tabulated by LF for this extended emission. From 
the contour levels and extent of the emission, our very rough 
estimate of the total flux density of N is 75 mJy, which we 
have indicated with a factor of 2 margin towards both lower 
and higher values in Fig. 7. 

Region N is also unmistakably present in the 0.61 GHz im- 
age of LF. Even though, for their own analysis, LF have in- 
cluded it as part of E, the third 0.61 GHz model component 
tabulated by LF (with flux density 0.11 Jy) clearly forms part 
of N, given its position relative to E. The model component 
might only represent part the total flux density from that area, 
and it is thus indicated as a lower limit in Fig. 7. 

At 0.33 GHz, LF have modeled the source as a single 
Gaussian, which they ascribed fully to E. However, in their im- 
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Fig. 7. Continuum component flux densities at multiple frequencies, 
based on data at 0.33, 0.61, 1.67, 2.3, and 8.4 GHz from LF, at 2.3, 
8.4, and 15.4 GHz from K03, and at 1.049 GHz from our own obser- 
vations. See Sec. 4 for details. 



age we find good evidence for the presence of all three com- 
ponents, W, E, and N. The image is clearly extended and elon- 
gated, and the western side is the fainter part. We believe this 
is W; all other data establish that W has a more inverted broad- 
band spectrum and is fainter than E at all frequencies below 
8.4 GHz (see Fig. 7). Unfortunately, E and N are difficult to 
disentangle in the 0.33 GHz image, both because of the lower 
VLBI resolution, and because of scatter-broadening, although 
LF may have overestimated the amount of interstellar scatter- 
ing at 0.33 GHz when they took the entire visible structure to 
be due to E. We find that the overall position angle of the con- 
tour image is 75°-85°. This suggests that the location of the 
centroid of the bright, eastern part of the source is influenced 
by significant contributions from both E (in PA~100° from W) 
and N (in PA~65° from W). Thus, we show in Fig. 7 a plau- 
sible but by no means unique decomposition of the total flux 
density of the single model component fitted by LF; to indicate 
our uncertainty, no marker symbols are drawn. 

Components W and E, we believe, are thus separately vis- 
ible at all frequencies, while component N is visible in the 
images up to 1.67 GHz. Above 1 GHz, where the linear res- 
olution is sufficient and the impact of interstellar scattering 
negligible, "component" E shows considerable substructure. 
At the highest observing frequencies, 8.4 GHz and 15.4 GHz, 
"component" W is seen to have compact substructure as well. 
Where measurable, the intrinsic angular sizes of W or its sub- 
components are considerably smaller than those of E or its sub- 
components. Also, from Fig. 7, W clearly has the most inverted 
low-frequency spectrum, and very likely also the highest peak 
frequency and the flattest high-frequency spectrum. 

N has quite the opposite properties: it is the most extended 
and does not appear to have much compact substructure; its 
broad-band radio spectrum is probably steep over the full fre- 
quency range (although with large uncertainties from the data 
of LF). Thus, W, E, and N, successively are in a sequence of 



Vermeulen et al.: Hi in the one-sided "compact double" radio galaxy 2050+364 



7 



decreasing compactness, steepening high frequency spectral in- 
dex, and decreasing peak frequency. 

5. Discussion 

5.1. 2050+364: a one-sided core-jet structure ? 

The three regions W, E, and N show a progressive change of 
compactness and of spectral shape which is not typical for com- 
pact double-lobed radio sources, but rather for core-jet radio 
sources, where relativistic beaming and opacity effects lead to 
a one-sided observed structure. Furthermore, the relative align- 
ment of the three regions is also compatible with a one-sided 
core-jet source. The very compact sub-components of W seen 
at the highest frequencies are aligned roughly, but not exactly, 
along the line towards E. The relative disposition of the sub- 
structures in E at various frequencies suggests that this may 
well be radio emission from shocks at the location of a sud- 
den bend in a single jet, which is first visible at W, and after 
bending at E continues towards the more extended emission re- 
gion N. The jet could be at a fairly small angle to the line of 
sight, such that the deprojected angle of the bend at E need not 
be large; this is seen in many other core-jet sources. Due to 
its substructure, E has a somewhat unusual convex broad-band 
radio spectrum, peaked at frequencies of a few GHz. We find 
that 2050+364 is set apart from many other GPS/CSS sources 
by the fact that at all GHz frequencies its total flux density is 
dominated by two emission regions, W and E, which are rather 
compact and widely separated, i.e., each one spans only a small 
fraction of the distance between them. Genuine double-lobed 
GPS/CSS sources can show compact hot spots at the two ends, 
but then in addition they typically have clearly visible, often 
dominant, extended emission regions (the lobes). Region N in 
2050+364, on the other hand, is rather faint, even at the lowest 
frequencies, and as such would be rather atypical if it were a 
lobe beyond its putative hot-spot E. A convincing counter-lobe 
is not seen. Classical core-jet sources, on the other hand, do of- 
ten show a faint extended emission region beyond a shock at a 
bend in the jet. 

We think that it is, furthermore, plausible that the sub- 
components of W may be at or adjacent to the traditional radio 
core. Many core-jet radio sources show a few bright knots in 
the innermost jet. The overall apparent separation velocity be- 
tween W and E at 1 .67 GHz has been delimited by LF to be 
no faster than c. However, this limit is still several times faster 
than what seems to be typical in CSOs (e.g., Gugliucci et al. 
2005), and it also does not rule out that some sub-components 
might even show superluminal motions, as is sometimes seen in 
other sources, where stationary radio components are thought 
to mark regions of either confinement or bending in a jet, while 
other radio features may be moving with a flow or due to trav- 
eling disturbances. It is possible that the sub-components of W 
are still part of the inner jet, with the true core being somewhat 
further to the west, perhaps at a more compact radio compo- 
nent, as yet unseen due to a spectrum which is significantly in- 
verted as a result of synchrotron self-absorption, free-free ab- 
sorption, or a combination of both, as is probably the case in 
several other compact radio sources (e.g., Kameno et al. 2003). 



It would be interesting to attempt to monitor the positions of 
some of the well-defined sub-components in the W and E re- 
gions at high frequency. 

In most well-known two-sided GPS/CSS sources, at the rel- 
atively high frequency of 15 GHz synchrotron self-absorption 
and/or free-free absorption do not fully hide the central core 
component (to a reasonable dynamic range of a hundred, say). 
However, even higher frequency observations would be useful 
before definitely ruling out the existence of a more inverted- 
spectrum core component between W and E, in which case the 
radio source might be two-sided after all. Likewise, even if the 
currently imaged structure is one-sided, with the core at or near 
W, it is possible that further high-dynamic-range imaging could 
reveal emission from a second jet or lobe on the opposite side. 
This situation has occurred, for example, in 2352+495, a well- 
known CSO (see Conway et al. 1992). 

The original selection of 2050+364 as a Compact Double 
(suggesting two-sided radio emission) seems to have been the 
result of a similarity between W and E which is partly coinci- 
dental; we now know that their properties diverge when studied 
over a wider range of frequencies, and that a section further out 
in the jet, N, has a steeper, more typical broad-band spectrum. 
Several other examples are now known in which the similarity 
between the components of an apparent compact double breaks 
down when the source is studied with high enough resolution 
and dynamic range and at high enough frequency. Amongst the 
sources now considered to be one-sided is CTD93 (see Shaffer, 
Kellermann, & Cornwell 1999); like 2050+364, it was a mem- 
ber of the original group of CDs. The identification of a cen- 
trally located core component is a crucial part of confirming 
whether a CD or GPS source is indeed two-sided, and there- 
fore a CSO, rather than a one-sided core-jet source in which 
two emission regions happen to have comparable flux densi- 
ties and sizes, and therefore roughly similar spectra (Conway 
et al. 1994, Wilkinson et al. 1994). 

5.2. The optical redshifts and H i: infall or outflow ? 

The optical spectroscopy reported by de Vries et al. (2000) 
gives z = 0.3547 consistently for both [Oin]/l4959 and 
[Om] /I5007, while H/3,14861 is at z = 0.3536. Thus, the ob- 
served centroid velocity of the [O m] doublet is formally about 
325 kms -1 larger than that of H/3. We believe that this differ- 
ence is probably significant. The spectral resolution was 15 A, 
but the centroid wavelengths are listed in integer A, i.e with 
4 significant figures. The [Om] doublet lines are narrow and 
clearly have an excellent signal-to-noise ratio, so the total un- 
certainty in the centroid of these lines could well be of order 
±1A in wavelength, or +0.0002 in redshift, or ±50 kms -1 in 
velocity. The observed FWHM of Hy6 looks like it is roughly 
50 A, or 2500 kms~', and the line has a lower signal-to-noise 
ratio, so the uncertainty in its centroid is probably several times 
larger than for the [O m] doublet. Nevertheless, we think that 
the combined velocity uncertainty probably cannot account for 
the entire 325 km s _1 observed centroid velocity offset between 
[Om] and H/3. 
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The [O m] line velocity and the velocity of the deepest H i 
features agree to within 10 kms~' (see Fig. 4). The radio spec- 
tral resolution is 4.5 km s _1 , and the H i absorption features are 
present at high signal-to-noise, so the uncertainty in this match 
is dominated by the optical data, and we believe that the offset 
of the [O m] velocity from that of the deep H i feature is not 
significant. 

The VLBI data show that there is atomic gas at the [O in] 
velocity spanning a projected distance of several hundred par- 
sees, and probably overlapping in projection with the radio 
core. This makes it likely that the [Om] emission is from a 
conventional narrow line region (NLR) surrounding the active 
nucleus, rather than from an isolated cloud moving with a fairly 
high peculiar velocity. In Sec. 5.3 some less conventional pos- 
sibilities are discussed, but we believe these are less plausible. 
A substantial centroid offset from the systemic velocity of the 
host galaxy seems in general more likely for the broad line re- 
gion (BLR), which shows the kinematics of gas under the im- 
mediate local influence of the active galactic nucleus, than for 
a region with an extent of several hundred parsecs. Thus, we 
have adopted the [O in] redshift, z = 0.3547, as the redshift of 
the host galaxy, and the zero point of the velocity scale used 
for our main analysis. This means that both the H/3 line and 
the lower opacity H i absorption features extend to negative ve- 
locities (i.e. outflowing, approaching the observer); in the rest- 
frame, the velocity offset of H/? is -244 km s _1 . 

The Hi absorbers in 2050+364 thus appear to be some- 
what different than in PKS 1549-79 (Tadhunter et al. 2001) and 
B1221-42 (Morganti, priv.comm.). In those objects the veloc- 
ity of the Hi line does not coincide with the optical redshift 
of the [O in] line, but rather with the [O n] line. Because of its 
lower ionisation level, one might expect [On] to originate in 
a larger region and perhaps to be more representative of the 
systemic velocity than [Om]. 

5.3. The locations of the H i absorbers relative to the 
radio source 

We think that all of the H i absorption may well occur in the 
same region as the radio source, in the inner kiloparsec of the 
host galaxy, and that the presence of this atomic gas might have 
some connection with the active nucleus, or with the radio mor- 
phology. Below we argue that an association of the narrow, 
deep Hi component with the NLR in 2050+364 is plausible, 
and that the broader H i absorption towards region W could be 
related to the BLR. 

The atomic gas is probably not kinematically disturbed 
very much as a result of interaction and bending of the jet at 
E. The absorption profile there only has two narrow compo- 
nents, and moreover, the dominant deep line is rather similar to 
the one at W, where no bending of the jet is visible. If the jet at 
E does impact and disturb clouds containing atomic hydrogen, 
these probably do not cover the radio source as seen along our 
line-of-sight. 

The sharp, deep absorption features at low velocity are at 
first sight the most puzzling to explain. The integrated peak 
opacity in 2050+364, 16 %, is the third deepest in the sample of 



41 CSS and GPS sources analysed by Pihlstrom et al. (2003). 
Based on the anti-correlation of peak opacity with linear size 
in that sample, < 1 % absorption would be more typical for 
a projected extent of ~300 pc. Conversely, the FWHM of the 
integrated deep H i absorption feature in 2050+364, 16 km s _1 , 
is the narrowest in the sample, where most of the FWHM are 
in the range of fifty to several hundred km s _1 , with no obvious 
dependence on source linear size. 

Because the remarkably high peak opacity is partly bal- 
anced by a remarkably low FWHM, the integrated column 
depth lies only mildly above the anti-correlation with linear 
size in the sample of Pihlstrom et al. (2003): assuming uni- 
form coverage and a spin temperature T sp ~ 100 K, as is of- 
ten thought to prevail in typical ISM conditions, we obtain 
N(Hi) ~ 5 x 10 20 (r s/ ,/100) cm 2 ; column depths for the in- 
dividual features are listed in Table 2. 

The most likely location for this H i absorption is in the 
neutral cores of the clouds in the NLR, given the excellent cor- 
respondence with the velocity centroid of the [O m] line, and 
the projected extent of at least 300 pc. 

Taking «hi = 100 cirT 3 as a rough estimate of the atomic 
gas density in NLR clouds, and using the 300 pc transverse 
extent covered by the absorbers as an estimate of the line-of- 
sight depth through the NLR as well, this would imply that 
clouds are present along about 1 % of the pathlength. The key 
distinguishing property, then, is the low kinematic dispersion. 
The main direction of motion and/or rotation in the NLR of 
2050+364 is evidently perpendicular to our line-of-sight, and 
coherent in velocity to ~10 kms~' over several hundred par- 
secs both along the line-of-sight and transverse to it. It is this 
velocity coherence which allows a deep, narrow absorption line 
to build up. Thus, it seems unlikely that the absorption is re- 
lated to directly inflowing gas, feeding the nucleus, because 
the velocity is so similar over a region substantially larger than 
typical accretion regions. 

On the other hand, the absorption features of comparatively 
low opacity but spanning a broader observed velocity range (to 
nearly -200 kms -1 ) occur only towards W, which we believe 
is most likely to be at or close to the nucleus (see Sec. 5.3). 
It is plausible that this shows neutral hydrogen which is either 
being entrained by the inner jet of 2050+364, or is flowing out 
from the accretion region. Perhaps this atomic gas is related 
to the BLR in 2050+364, particularly since the H/3 centroid, 
while subject to considerable uncertainty, is probably also at 
a negative velocity of a few hundred kms -1 . The neutral gas, 
especially that part of it which happens to be visible in ab- 
sorption, is likely to sample only a part of the velocity profile 
of the ionised gas. The broader part of the Hi profile clearly 
has several distinct kinematic components (a possible repre- 
sentation with Gaussian lines is given in Fig. 6 and Table 2); 
the average opacity is t ~ 0.025. Assuming W to be uni- 
formly covered gives a column depth integrated over the ob- 
served velocity range -50 kms" 1 to -180 kms -1 , of roughly 
N(H I) ~ 6x l0 2Q (T sp / 100) cm 2 . Making the further and surely 
oversimplified assumption that the absorbing atomic gas is also 
uniformly dense along a pathlength comparable to the total 
transverse extent of W (~10 pc), then the atomic gas density 
would be n H i = 20 cirT 3 , which is rather low. Unfortunately, 
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with the EVN we cannot resolve region W at the frequency of 
the H i line. The optical BLR itself probably has an extent of at 
most 1 pc, so if the atomic gas were also confined to that region 
only, then its density would be perhaps two or three orders of 
magnitude higher, because the pathlength would decrease, the 
opacity and column depth would increase since the covering 
factor of W would be less than unity, and finally, close to the 
active nucleus the spin temperature could be significantly ele- 
vated. 

Could the Hi absorption arise instead outside the host 
galaxy of 2050+364 ? A chance superposition, in which the 
galaxy visible at the position of 2050+364 would be a fore- 
ground object rather than the host of the radio source, is a priori 
unlikely. Such a situation does exist for lensed radio sources, 
such as 0218+357, where the foreground lensing galaxy also 
leads to Hi absorption (Carilli, Rupen, & Yanni 1993), but we 
do not think that the radio morphology of 2050+364 is sugges- 
tive of gravitational lensing. 

Close to the radio position, two other galaxies of compara- 
ble brightness are visible in the images of Biretta et al. (1985 
and O'Dea et al. (1990). One or both of these could be a true 
companion to the host galaxy of the radio source. The closer 
neighbour, in PA 180°, is centred at a projected distance of only 
2 arcsec (10 kpc at z = 0.3547), and so part of its disk could 
well be overlapping with the radio source in projection. The 
other one, centred at a projected distance of 7 arcsec (35 kpc at 
z = 0.3547) in PA -120°, seems to show an extension (a spiral 
arm ?) in the direction of the radio source, and this might be 
another candidate foreground absorber. 

However, in any external absorber scenario considered, that 
other galaxy would probably have to be responsible for both 
the [O in] line as well as the H i absorption, given the velocity 
correspondence to within 10 kms~'. But that would mean that 
the other galaxy may also have an active nucleus, or prodigious 
star-bursting activity, since strong [O m] emission is not usually 
seen in ordinary galaxies. The optical spectrum shown in de 
Vries et al. (2000) suggests a high ratio of [Oin] to H/3, i.e., 
high ionisation gas. This would indicate ionisation by an AGN 
rather than a starburst. In our view the required coincidence 
of two active galaxies makes external absorber options much 
less plausible than the simple model discussed earlier that the 
Hi absorption shows atomic hydrogen in the NLR and BLR 
associated with the compact radio source 2050+364. 

6. Summary 

We have presented and discussed VLBI continuum and spectral 
line data at 1049 MHz of the compact radio source 2050+364, 
and we have interpreted the multi-frequency continuum VLBI 
images of LF. Our conclusions can be summarised as follows: 

- The continuum structure at 1049 MHz consists of a faint 
component N to the northeast, in addition to the two well- 
known components E and W, which constitute the original dou- 
ble; their flux density ratio is 3:1. 

- The compact-double properties hold only over a limited 
frequency range. W, E, and N, in that order, have decreasing 
compactness, steepening low frequency spectral index, and de- 
creasing peak frequency. 



- Compactness and radio spectra, plus the alignment of 
substructures within W and E, suggest that 2050+364 is a one- 
sided core-jet originating at or near W, bending sharply (in pro- 
jection) at E, and continuing towards N. 

- At the high velocity end of the absorption spectrum is a 
deep H i line, which reaches a maximum opacity t = 0.26 at 
E, with a FWHM of 13 km s _1 . At W it is even narrower, with 
a FWHM of 10 kms -1 and somewhat lower opacity, r = 0.21. 
At N the opacity is t ~ 0.09. This absorption thus covers the 
entire source, extending over > 300 pc, albeit with significant 
differences in opacity. Implied column depths are 7V(H i) ~ 5 x 
10 20 (r sp /100)cm- 2 . 

- Extending to lower velocities is lower opacity Hi ab- 
sorption, which only covers W. The average opacity between 
-50 kms -1 and -180 kms -1 is r ~ 0.025, and the estimated 
column depth is 7V(H I) ~ 6 X 10 20 (r s/J /100) cm" 2 . 

- The centroid of the [Om] optical doublet lines, z = 
0.3547, coincides to within 10 kms~' with the distinct deep 
line at the high velocity end of the Hi radio absorption line 
spectrum. The formal line centroid of the optical Hfi line is at 
an observed velocity of -325 km s _1 relative to [Om]. 

- We believe that the uncommonly deep but also uncom- 
monly narrow H i absorption line is likely to be due to atomic 
gas in the cores of NLR clouds in the inner kpc of 2050+364. 
Assuming NLR clouds of density «hi = 100 cm -3 in a region 
with a 300 pc radius would imply their presence along 1 % of 
the line-of-sight. The direction of motion and/or rotation in the 
NLR of 2050+364 is coherent to ~ 10 km s _1 over several hun- 
dred parsecs both along the line-of-sight and transverse to it, 
and is largely perpendicular to our line-of-sight. 

- We believe that the lower opacity absorption at W ranging 
to -200 km s _1 is likely to be due to atomic gas which is either 
being entrained in the inner parsecs of the jet of 2050+364, or 
is flowing out from the accretion region. It is plausibly related 
to the BLR in 2050+364. The atomic gas density is at least 
«Hi = 20 cm -3 , for a 10 pc absorbing region, but could well 
be as much as three orders of magnitude higher if the region is 
smaller. 
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